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PONDEROMOTIVE  EFFECTS  IN  NONNEUTRAL  PLASMAS 


B.  M.  Lamb  and  G.  J.  Morales 
Physics  Department 
University  of  California 
Los  Angeles,  California  90024 


ABSTRACT 


This  investigation  considers  the  ponderomotive  effects  which  arise  in  a 
single  species  plasma,  i.e.,  a  nonneutral  plasma.  The  important  difference 
from  a  neutral  plasma  is  that  quasi-neutral  density  cavities  given  by 
6n/nQ  =  —  I E [  /16irnoT  cannot  arise  in  a  single  species  plasma.  Instead,  it 
is  found  that  the  ponderomotive  force  is  balanced  by  self-consistent  space 
charge  fields,  and  results  in  Sn  »  V  |E„|  /16irmii>  ,  where  u>  is  the  fre¬ 
quency  and  E„  refers  to  the  wave  electric  field  parallel  to  the  confining  mag¬ 
netic  field.  In  addition  to  the  density  rearrangement,  the  zero  order  E  *  B 
rotation  is  modified.  The  self-consistent  nonlinear  state  has  been  found  for 
a  large  amplitude  £  =  0  axially  standing  Gould-Trivelpiece  wave.  The  linear  dis 
persion  relation  is  modified  by  the  presence  of  the  large  amplitude  wave,  and 
results  in  a  nonlinear  frequency  shift.  This  shift  produces  an  interesting 
hysteresis  effect  in  the  nonlinear  resonant  response  of  the  plasma  when  the 
frequency  of  the  external  driver  is  swept  slowly. 
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I .  INTRODUCTION 


Nonlinear  wave  effects  in  plasmas  fall  into  two  broad  categories:  wave- 
particle  interactions  and  fluid  nonlinearities.  For  waves  whose  phase  velo¬ 
cities  are  much  larger  than  the  average  thermal  velocities  of  all  particle 
species,  the  wave-particle  interaction  plays  an  insignificant  role.  In  such 
systems  the  fluid  nonlinearities  are  primarily  responsible  for  the  modifica¬ 
tion  of  the  zero  order  plasma  equilibrium  as  the  wave  amplitude  is  increased 

''  /  ' 

beyond  the  linear  regime.  The  subject  of  the  present  analytical  study  falls 
into  the  latter  class. 


Among  the  numerous  fluid  nonlinearities  associated  with  a  large  amplitude, 
high  frequency,  electrostatic  wave,  there  are  two  effects  of  overwhelming  impor¬ 
tance:  parametric  instabilities  and  ponderomotive  force  density  modifications. 

In  conventional  neutral  plasmas  both  of  these  can  appear  simultaneously  and 
lead  to  interesting  physical  effects  which  have  been  extensively  investigated 
theoretically  and  experimentally.  However,  in  single  species  plasmas,  i.e., 
nonneutral  plasmas,  the  parametric  instability  channel  is  not  readily  available 
because  there  is  only  a  single  mass  present.  Consequently,  in  these  systems 
it  is  expected  that  the  first  nonlinear  effect  to  appear  is  modification  of 
the  zero  order  density  by  the  ponderomotive  force.  This  is  the  specific  issue 

that  is  considered  in  this  work. 

1  2 

It  is  widely  accepted  *  that  in  neutral  plasmas  the  density  changes  6n 

produced  by  the  ponderomotive  force  are  given  by  the  expression 

Sn/n  =  —  I E I  /16irn  T,  where  n  is  the  zero  order  density,  T  the  plasma  tempera- 
o  1  ^ 1  o  o 

ture,  and  |e|  the  amplitude  of  the  high  frequency  field.  Physically,  the  end 
result  can  be  viewed  as  a  pressure  balance  between  the  wave  and  the  plasma. 
However,  in  the  calculation  leading  to  such  a  result,  the  existence  of  a  neu¬ 
tralizing  species  is  essential.  The  sequence  of  events  can  be  understood  by 


the  following  time  ordering.  First  the  high  frequency  field  acts  on  the  lighter 
species  and  pushes  these  particles  out  of  the  regions  where  | e(  is  maximum.  An 
intermediate  space  charge  separation  results  from  the  relative  displacement  of 
the  two  species.  Finally,  the  heavier  species  responds  to  the  space  charge  field 
and  follows  the  lighter  species,  thus  resulting  in  a  density  depression  which  is 
overall  neutral.  From  this  picture  it  is  apparent  that  when  the  second  species 
is  absent,  the  resulting  density  change  is  significantly  modified.  In  particular, 
balancing  kinetic  pressure  against  field  pressure  is  no  longer  possible.  It  is 
the  aim  of  this  study  to  provide  a  relevant  description  of  the  ponderomotive 
effect  when  the  neutralization  due  to  a  second  charge  species  is  absent. 

Technically,  the  difficulty  in  calculating  the  ponderomotive  effect  in  a 

nonneutral  plasma  arises  because  of  the  need  to  satisfy  a  global  self-consistent 

equilibrium  in  finite  geometry.  In  a  nonneutral  system  once  a  charge  separation 

is  induced,  all  zero  order  quantities,  e.g. ,  density,  space  potential,  and  rotation 

are  modified.  In  addition,  the  strong  symmetry  requirements  imposed  by  the  zero 

order  confinement  scheme  must  be  incorporated.  To  obtain  a  concrete  result  that 

includes  all  these  effects,  the  present  study  isolates  the  ponderomotive  changes 

3  U 

produced  by  an  axially  standing  Gould-Tr ivelpiece  mode  *  in  a  nonneutral 
plasma  column  that  is  strongly  magnetized.  Aside  from- defining  the  prototype 
problem  in  this  class,  the  situation  is  realizable  in  the  various  nonneutral 
plasma  devices  developed  recently^’^  which  exhibit  excellent  confinement  proper¬ 
ties. 

The  principal  result  of  this  study  is  that  a  fully  self-consistent  pondero¬ 
motive  change  can  be  obtained  analytically  for  the  case  of  an  azimuthally  sym¬ 
metric  (9  =  0)  Gould-Trivelpiece  mode.  Near  the  center  of  the  column  the  den¬ 
sity  change  is  found  to  be  proportional  to  the  Laplacian  of  the  field  amplitude, 


2  2  2 
i.e.,  fin  a  V  |  E|  instead  of  the  conventional  result  6n  ®  -|  E|  for  a  neutral 

plasma.  In  addition,  the  consequences  of  the  density  change  have  been  calculated. 
It  produces  a  nonlinear  frequency  shift  of  all  modes  satisfying  the  proper  wave 
number  matching  conditions.  An  interesting  hysteresis  effect  is  found  in  the 
resonant  response  of  the  standing  Gould-Trivelpiece  mode.  It  arises  when  the 
nonlinear  frequency  shift  is  comparable  in  magnitude  to  the  collisional  broad¬ 
ening  of  the  mode.  A  similar  effect  has  been  obtained  in  a  recent  low  tempera¬ 
ture  experiment^  involving  a  two  dimensional  nonneutral  plasma.  It  is  found  here 
that  a  three  dimensional  nonneutral  plasma  can  also  exhibit  this  nonlinear  phe¬ 
nomenon. 

The  manuscript  is  organized  as  follows:  In  Sec.  II  we  derive  the  pondero- 
motive  effects  due  to  a  £  =  0  axially  standing  Gould-Trivelpiece  wave  without 
specializing  to  a  specific  zero  order  equilibrium.  In  order  to  obtain  concrete 
results  we  consider  a  rigid  rotor  equilibrium  in  Sec.  III.  In  Sec.  IV  the  non¬ 
linear  frequency  shift  produced  by  the  ponderomotive  force  is  presented. 
Conclusions  are  given  in  Sec.  V. 


T.  DERIVATION  OF  PONDEROMOTIVE  EFFECTS 

Consider  a  cvlindricaliy  symmetric  nonneutral  plasma  consisting  of  particles 
of  charge  q  and  mass  m  confined  inside  a  conducting  cylinder  of  radius  a  by  a 
strong,  uniform  axial  magnetic  field  B  =  Boz,  as  illustrated  in  Fig.  1(a).  Cylin¬ 
drical  coordinates  are  utilized,  with  the  z  axis  coinciding  with  the  cylinder  axis. 
For  simplicity,  the  plasma  is  assumed  to  be  infinite  in  axial  extent.  This  assump¬ 
tion  removes  some  interesting,  but  complicating,  physical  effects  associated  with 
the  nonlinear  lengthening  which  can  occur  in  a  finite  length  plasma  column.  An 
axisymmetric  (£  =  0)  Gould-Trivelpiece  standing  wave  is  envisioned  to  be  resonantly 
driven  so  that  a  steady  state  amplitude  results  due  to  weak  collisions  with  neutrals 
The  problem  under  investigation  in  this  study  is  to  find  the  nonlinear  modifi¬ 
cations  of  the  zero  order  equilibrium  quantities,  e.g.,  density,  space  potential, 
and  rotation,  caused  by  a  large  amplitude  standing  wave.  The  effect  of  these 
changes  in  the  equilibrium  quantities  upon  the  properties  of  the  Gould-Trivelpiece 
modes  is  also  investigated.  Since  the  problem  is  periodic  in  the  axial  direction, 
it  need  only  be  analyzed  within  one  wavelength  of  the  standing  wave. 

The  essential  physics  of  the  problem  can  be  treated  using  a  fluid  description 
provided  that  the  phase  velocity  is  large  compared  to  the  thermal  velocity.  For 
Gould-Trivelpiece  modes  with  to  <  <<  ft  this  condition  is  readily  met.  Unlike 

theoretical  treatments  of  neutral  plasmas  which  are  often  modeled  as  being  infinite 
and  uniform,  nonneutral  plasmas  are  necessarily  of  finite  radial  extent  because  of 
the  resulting  large  self-electric  fields.  Due  to  this  property  one  must  consider 
a  specific  equilibrium  configuration  in  order  to  pose  a  well  defined  problem; 
however,  in  spite  of  this  restriction,  the  results  obtained  in  the  present  study 
are  of  a  fairly  general  nature. 

The  response  of  the  system  is  described  by  the  equation  of  motion, 
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'  jW  *y-  m '*  +  n/y* *  -  s  ^ 


the  continuity  equation, 


~ct  +  v  •  ci  ]0  =  0, 


and  Poisson's  equation, 


V  $  =  -4irq/}, 


where  =  qBQ/mc  is  the  cyclotron  frequency  v'V(r,  z,  t)  is  the  total  plasma 
density ,  ^(r,  z,  t)  is  the  fluid  velocity, ^(r,  z,  t)  is  the  thermal  pressure,  and 


'Kr,  z,  t)  is  the  potential.  The  self-magnetic  field  generated  by  the  plasma  rota¬ 


tion  has  been  neglected  in  Eq.  (1)  because,  for  the  plasma  parameters  considered, 


i.e.,  u)p  a/f2c  <<  1  (where  oj^  is  the  plasma  frequency),  it  is  small  compared 


to  Bq.  We  also  assume  the  electromagnetic  wavelength  is  long  compared  to  the 


characteristic  dimensions  of  the  physical  system,  i.e.,  u>a/c  <<  1,  so  that  the 


quasi-electrostatic  approximation  is  valid. 


We  proceed  to  solve  Eqs.  (1)  -  (3)  by  ordering  according  to  the  amplitude 


of  the  Gould-Trivelpiece  wave.  The  linear  wave  quantities  are  considered  to  be 


a  first  order  perturbation  to  the  zero  order  equilibrium  quantities.  In  second 


order  the  ponderomotive  force  appears  through  the  time  independent  (dc)  component 


f  i  '  I 

of  the  t  •  V  }  term.  At  this  stage  we  consider  only  the  dc  second  order  modi¬ 


fications  to  the  zero  order  equilibrium  quantities.  In  Sec.  IV  we  feedback  the 


second  order  modified  equilibrium  quantities  into  a  nonlinear  dispersion  relation 


to  obtain  a  nonlinear  frequency  shift  of  the  standing  wave.  The  effect  of  second 


order  contributions  at  harmonics  of  the  linear  wave  frequency  are  not  considered 


in  this  manuscript. 


The  electric  potential,  density,  velocity,  and  pressure  are  separated  into 


three  terms  consisting  of  zero  order  equilibrium,  linear  response,  and  dc  quadratic 


>  f 


a  4 

••  v  •  ‘i 


9  9i 


v-vv-1 


nonlinear  changes  nroduced  by  the  large  amplitude  wave: 


$  =  4>o  +  2"!$  exp(-ioit)  +  c.c.]  +  5<)>, 


i  .•»  i  - 

c'}  =  no  +  -jl11  exp(-ioit)  +  c.c.]  +  fin. 


•  •  1 

,  =  V  0  +  ^[v  exp(-ia)t)  +  c.c.]  +  5v, 

o  Z  _ 


Vj^  =  TV  (n  +  fin) , 


where  <|>o(r),  nQ(r),  and  V  (r)  are  zero  order  equilibrium  quantities,  4>  (r,z), 
n(r,z),  and  v(r,z)  are  complex  amplitudes  representing  the  linear  response  to 
the  high  frequency  wave  at  frequency  a),  and  6<j>(r,z),  <Sn(r,z),  and  <5v(r,z)  are 
the  steady  state  self-consistent  nonlinear  changes  produced  by  the  wave.  Constant 
temperature  is  used  for  the  dc  quantities  because,  as  we  later  find  in  Sec.  Ill, 
in  the  bulk  of  the  plasma  column  temperature  plays  a  small  role  in  the  force 
balance  against  the  ponderomotive  force.  In  addition,  no  explicit  usage  of  a 
collision  frequency  is  made  at  this  stage;  the  wave  amplitude  is  simply  assumed 
to  be  constant.  Temperature  corrections  to  the  wave  quantities  are  neglected 
because  for  Gould-Trivelpiece  waves,  both  the  phase  velocity  and  the  group  velocity 
are  much  greater  than  the  thermal  velocity. 

An  equation  for  the  zero  order  equilibrium  is  obtained  by  substituting  from 
Eq.  (4)  into  Eqs.  (1)  -  (3)  and  setting  the  wave  amplitude  to  zero.  For  this  con¬ 
dition  the  r  component  of  Eq.  (1)  becomes 


iv2 

r  o 


-  -3-  <})  +  Q  V  -  —  ln(n  /n  )  , 

m  3r  o  o  m  3r  op 


in  which  the  density  has  been  normalized  to  the  density  on  axis  n^  =  n0^r  =  • 

1  3 

Operating  on  Eq.  (5)  with  —  —  r  and  using  Poisson's  equation  to  eliminate 


yields  the  differential  equation 


Xn27-|rrl—  ln(n/n)-n/n  =  — ^  (V  (Sir  +  V  )  ]  , 

D  r  9  r  3r  op  op  to  2  r  9  r  o  o 

P 

2  2  2  2  2 
where  w  =  4nn  q  /m  and  XQ  =  1/Kp  =  T/4nn  q  .  Using  the  scaling 


(6) 


p  5  r/XD,  N (p )  =  nQ(r/XD)/np,  A(p)  = 


V  (r/X  ) 
o _ D 

/T7m 


(7) 


gives 


1  9__ 
P  3p 


3p 


J.n  N 


-»»  -  r  [ A(flp / (o  +  A)] 
p  9p  p 


(8) 


The  boundary  conditions  for  Eq.  (8)  are  N(p  =  0)  =  1  and  N (p  =  aAp)  =  0. 
Equation  (8)  may  be  satisfied  by  assuming  a  given  form  for  the  density  profile 
N(p)  and  solving  for  the  resulting  scaled  velocity  A(p),  or  conversely  as  is 
done  in  Sec.  Ill,  by  taking  a  rigid  rotor  velocity  profile  and  solving  for  N(p). 

The  requirement  that  N  must  vanish  at  the  wall  places  constraints  upon  the  rota- 

/ 

tion  velocities  which  are  physically  allowed. 

The  linear  response  quantities  <}>(r,z),  n(r,z),  and  v(r,z)  are  separated 
from  Eqs.  (l)-(3)  giving  the  equation  of  motion, 

-iw  v  -  -  V  v.r+iv  4—  (rV  )0  =  -  ^  7  if  +  fi  v  x  z,  (9) 

_roo  r  r  9r  o  m 

A  A 

where  r  and  0  are  unit  vectors,  the  continuity  equation, 


-i(d  n  +  7  •  (n  v)  =  0, 
o  „ 

and  Poisson's  equation, 

V2  <f>  =  -  4irqn. 

Equation  (9)  is  solved  for  v  in  terms  of  <f>  giving 


(10) 


(11) 
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1C0  „ 

V  =  -  X  *  V 

-  qno  ; 


(12) 


in  which  the  susceptibility  x  is  defined  in  terms  of  the  dielectric  tensor 
through  e  =  1  +  4ttx*  The  components  of  the  dielectric  tensor  are 


e  =  I  e 


rr 

er0 

°  \ 

0r 

£O0 

°  ) 

(13) 

i 

0 

e  / 

ZZ  , 


in  which 


£rr  £00 


er0  “  £0r 


1  - 


“p  (r) 


w  -  [fi  +  2V  (r)/r](n  +  1/r  3/3r(rV  (r))] 
o  o 


-i  [0  +  1/r  3/3r  (rV  (r))]w  (r) /to 

_ 2 _ E _ 

to2  -  [n  +  2V  (r)/r][n  +  1/r  3/3r(rV  (r))] 
o  o 


e  =  1  - 

zz 


to  2  (r) 
P 

2 

to 


(14) 


Equation  (12)  is  combined  with  Eqs.  (10)  and  (11)  giving  a  partial  differential 
equation  for  4> 


V*  e*V<t>=0.  (15) 

The  linear  normal  modes  of  the  nonneutral  plasma  system  are  obtained  by  solving 
Eq.  (15)  subject  to  the  boundary  condition  that  <J>(r=a,z)  =  0  and  the  condition 
that  <(>  is  finite  on  axis. 

The  time  averaged  (dc)  nonlinear  quantities  which  are  quadratic  in  the  wave 
amplitude  are  separated  from  Eqs.  (1)  and  (3)  giving 

-  —  V  SVn  r  +  q  V  <5<J>  -  mTJ  <5V  r  +  T  V  6n/n  =  -  m/2  Re(v*  •  V  v) ,  (16) 

ro0n  0  o  _ 
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E.- 

k'' 


e- 

a 


t.-'- 


ts 
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*.  "  \  •.  7*  7*  /'  .,*•  ,*• 


and 


V  64>  =  -4irq  5n. 


(17) 


Notice  that  Eq.  (16)  contains  the  driving  term  -m/2  Re(v*  •  V  v) ,  which  is  the 
ponderomotive  force  of  the  wave.  The  factor  of  one-half  results  from  taking 
the  time  average.  Substituting  Eq.  (12)  into  -m/2  Re(v*  ■  V  v)  gives 


-m/2  Re(v*  •  V  v)  =  - 


.222 
4tt  q  oi 


j- 

*  2 

3 

l  Z,  |  2 

9  i;  1 2 

;  z 

X 

3  z 

E  +  y  y 

•z—  E 

zz 

,  z .  rr  zz  dz  1  r  j 

(18) 
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in  which  E  and  E  have  been  defined  using 
r  z 


E  r  +  E  z  =  -Vi}>, 
r  z 


(19) 


and  derivatives  of  x/wp  have  been  neglected  since  they  are  of  higher  order 

2  2 

in  Up  /ft  ,  which  is  assumed  throughout  this  work  to  be  a  small  correction.  Since 


we  are  considering  a  strongly  magnetized  plasma  with  u  <  <<  ft,  the  general 


2,„2 


ponderomotive  force  in  Eq.  (18)  is  approximated  as  the  lowest  order  term  in  u^  /ft  , 

2 


-  m/2  Re (v*  •  V  v)  -  - 


3  *  lEz' 


3  z  ,  2 

4mw 


(20) 


which  physically  means  that  only  the  component  of  the  ponderomotive  force 
parallel  to  the  magnetic  field  is  retained.  With  this  approximation  the 


z  component  of  Eq.  (16)  becomes 
6<J>  + 


3_ 

3z 


a.  6*  +  q +  1  An 

m  9  .22  m  n 

4m  u  o 


=  0, 


(21) 


which  states  that,  in  the  strongly  magnetized  limit,  the  self-consistent  elec¬ 
trostatic  force,  modified  pressure  force,  and  ponderomotive  force  must  be  bal¬ 
anced  along  each  magnetic  field  line.  This  equation  implies  that  the  quantity 
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-  J 


A 

I  :tm 
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•  .  -  •  N 


Lfc.1.  mJL  M 


.  1 
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within  the  brackets  may  be  set  equal  to  a  function  f(r),  independent  of  z, 

.2, ~  ,2 

T 

(22) 


q  q  lEzl  T 

m  +  TTT  +  m  6n/no  =  f(r)< 


4ra  u 


Any  f(r)  which  vanishes  at  r  =  a  represents  an  allowable  equilibrium.  However 
in  the  context  of  a  deterministic  experiment,  one  must  find  the  appropriate  f(r) 
which  corresponds  to  the  equilibrium  which  the  system  evolves  to  as  the 
wave  amplitude  is  increased  adiabatically  in  a  well  defined  zero  order  non¬ 
neutral  plasma.  Two  important  constraints  are  used  in  determining  f(r). 

Particle  number  should  be  conserved  because  an  l  =  0  wave  in  a  strongly  magne¬ 
tized  plasma  should  not  cause  radial  plasma  transport,  and  because  in  our  axi¬ 
ally  infinite  model  particles  may  not  be  lost  axially.  Of  course,  in  an  experi¬ 
ment  particles  are  axially  confined  by  large  electrostatic  barriers. The 
second  constraint  is  that  canonical  angular  momentum  must  be  conserved  because 
the  excitation  of  an  H  =  0  wave  does  not  apply  an  external  torque  to  the  system. 
The  statement  of  particle  number  conservation  is 


\ 

2 


r+X/4 

dz 


-X/4 


•a 

rdr  6n  =  0, 
o 


(23) 


2,„2 


and  the  conservation  of  canonical  angular  momentum  P0  to  order  0)^“/^ 


is 


f+X/4 

■a 

1  \ 

■+X  /  4 

ra  2 

dz 

2irrdr 

qrAQ6n/c 

=  2 

dz 

2iTrdr  r  6n  =  0. 

(24) 

-X/4 

o 

4 

-X/4 

o 

The  z  integrations  in  Eqs.  (23)  and  (24)  are  performed  over  one-half  wavelength  X 
of  the  axial  standing  wave  because  6n  is  expected  to  have  axial  variations 
with  period  X / 2 .  Both  the  above  conditions  can  be  satisfied  if 


X_ 

2 


•+X/4 

6n  dz  =  0 
-X/4 


(23) 
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which  means  physically  that  plasma  is  rearranged  along  each  field  line  in  re¬ 


sponse  to  the  ponderomotive  force  of  the  standing  wave  without  transport  across 


field  lines. 


We  obtain  f(r)  by  operating  with  V  on  Eq.  (22),  substituting  Poisson's 


equation. 


-  V2  6n/n  -  ^3-  6n  - - 3_^  V2  |e  |  2  +  V2  f  (r)  , 

m  om  ,22  1  z1 

4m  to 


and  integrating  in  z  from  -X/4  to  +X/4: 


£(r)  ■  rir I  iE/  dz- 

4"“  J-X/4 


in  which  use  has  been  made  of  the  sinusoidal  axial  variation  of  6n  and  E 

z 


Having  determined  the  correct  choice  for  f(r),  Eq.  (22)  is  combined  with 


Eq.  (17),  giving  a  partial  differential  equation  for  6i}>, 


X,2  V2  =  -3-^  6|Ej2, 

4mto 


in  which 


I*/ -I 


and  having  the  boundary  condition  5<f>  (r  =  a,  z)  =  0.  One  may  interpret  Eq.  (28) 


by  saying  that  the  term  -I/AttX^2  (q /4rao  )  6  [  E^  |  represents  an  effective  nonlinear 


charge  density  induced  by  the  wave  which  is  then  Debye  shielded  by  the  plasma. 


In  the  appropriate  plasma  limit  for  which  KQa  >>  1,  the  V  6<t>  term  of 


Eq.  (28)  becomes  small  compared  to  6<J>  and  one  may  approximate  6<J>  as 


6<t>  a  - 


T  6'Ezl 


t- 
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Substituting  Eq.  (27)  into  Eq.  (26)  results  in  the  following  partial  differ¬ 
ential  equation  for  6n/n 

o 


A  2  V2  5n/n  -  N(r/An)  6n/n  - - — r  V2  6 1 E  |2.  (31) 

U  O  L)  O  ,  ,  Z  Z 

16irn  into 

p 


The  correct  boundary  condition  for  Eq.  (31)  when  T  ^  0  is  seen  from  Eq.  (21)  to  be 

6n/nQ  ->  0  as  r  a  because  both  6<j>  and  E^  are  zero  on  the  boundary. 

In  the  bulk  of  the  plasma,  i.e.,  outside  of  a  sheath  region  near  r  =  a, 

2  2 

the  first  term  of  Eq.  (31)  goes  as  A^  /a  ;  in  this  region  <5n  may  be  approximated 
as 

<5n  - - — j  V2  5  |  E  |2.  (32) 

16irmo)  Z 


Finally,  the  azimuthal  nonlinear  velocity  modification  6VQ  is  obtained 
from  the  radial  component  of  Eq.  (16), 


6ve 


1  3_ 
3r 


(33) 


in  which  we  have  used  the  definition  n*  =  n  +  2VQ(r)/r  of  the  effective  cyclo¬ 
tron  frequency  in  the  rotating  plasma  frame. 

Substituting  Eqs.  (22)  and  (27)  into  Eq.  (33)  results  in 


<5V 


0 


4  m  a)  ft* 


(34) 


We  note  that  6V  does  not  depend  explicitly  on  the  plasma  temperature  T, 
even  though  it  is  the  sum  of  the  E  x  B  drift  resulting  from  —  q  V  and  the 
diamagnetic  drift  resulting  from  -  T  V  <5n  with  the  magnetic  field  being  the 
effective  magnetic  field  in  the  rotating  plasma  frame.  The  physical  interpre¬ 
tation  of  Eq.  (21)  is  that  the  electric  force  and  the  pressure  force  combine  to 


ITT.  Rigid  Rotor  Equilibrium 


Q 

Most  treatments  of  nonneutral  plasma  equilibria  specialize  to  rigid  rotors,  * 

i.e.,  V  (r)  =  co  r  such  that  the  plasma  executes  a  bulk  rotation  about  its  axis. 

O  K 

In  this  section  we  consider  a  rigid  rotor  with  a  density  profile  which  is  flat 
from  r  =  0  to  within  several  Debye  lengths  of  r  =  a.  Equation  (8)  reduces  to 

^  fe-  °  If  N  -n  +  1-y-  (35) 


with 


2cor(<or  +  n) 

Y  9 - (o~2 - L- 

P 


(36) 


This  equation  has  been  previously  solved.  For  y  -  0  it  possesses  solutions 

for  which  N  is  uniform  for  many  Debye  radii  and  smoothly  falls  to  zero  in  a  sheath 
region  of  several  Debye  lengths  in  radial  extent,  as  shown  in  Fig.  1(b)  for  y  =  .0 
and  a  =  10 

The  approximate  values  of  u)„  allowed  by  Eq.  (35)  are  found  by  setting  y  =  0 

in  Eq.  (36)  and  solving  for  w  .  This  gives 

K 


U-  -  -  %  1  ±  /l  -  2u  2/n2 

K  Z  p 


(37) 


The  slow  rotational  mode  given  by  the  minus  sign  in  (37)  ,  in  the  strongly  magne- 
2  2  2 

tized  limit  of  u  <<  0  ,  is  (»)„  -  -w  /2f!,  which  is  simply  the  E  *  B  velocity 

p  K  p  «... 

due  to  the  associated  space  charge  field. 

The  solutions  of  Eq.  (15)  for  w  <  u  <<  Q  are  the  Gould-Trivelpiece  modes 

P 

of  a  nonneutral  plasma  column.  For  the  purpose  of  studying  these  waves  we  shall 
neglect  the  sheath  region  by  approximating  the  zero  order  density  as  constant. 
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Since  <p  must  vanish  at  r  =  a,  this  is  a  good  approximation  if  the  sheath  region 


is  small  compared  to  the  column  radius  a.  A  constant  density  solution  of  Eq.  (15) 


with  <f>  =  0  on  the  conducting  boundary  at  r  =  a  is 


=  E  /k  J  (x  r/a)  sinkz, 
o  o  on  ’ 


in  which  E  is  an  electric  field  magnitude,  J  is  a  Bessel  function,  and  x  is 
o  o  on 


the  n  zero  of  the  J  Bessel  function.  The  dispersion  relation  becomes 

o 


.2  2 

k  a  =  - 


.  2.  2 

1-0)  /O) 

P 


If  the  phase  velocity  w/k  is  large  compared  to  the  thermal  velocity  /f/'m 


2  2  2  2 

or  equivalently  if  k  <<  o)  ,  Landau  damping  is  negligible.  Using  Eq.  (39) 


2  2  2  2  2  2  2 

this  requires  that  w  /w  <<  1  -  x  /K„a  •  Recalling  that  K_a  >>  1  for  typical 

p  on  I)  D 


nonneutral  plasmas,  we  see  that  in  the  regime  we  are  considering,  w  <  <<  Q, 


weakly  damped  waves  can  exist.  Thus  by  externally  driving  the  plasma  at  frequency 


o),  the  excited  Gould-Trivelpiece  waves  can  become  steady  state  axially  standing 


waves . 


Using  the  defining  relation  for  6  Ez|  ,  Eq.  (29),  gives 


6 | E  | ^  =  —  E^  (x  r/a)  cos2kz. 
1  z1  2  o  o  on. 


For  the  first  radial  eigenmode,  i.e.,  n  =  1,  this  quantity  decreases  in  the 


radial  direction  as  shown  in  Fig.  2,  and  varies  axially  as  cos2kz. 


Inserting  Eq.  (40)  into  Eq.  (31)  gives 


V^(fin/n  )  -  (6n/n  )  =  -  ^5—  V^fj  ^(x  r/a)  cos2kz  -  (1  -  n  fn  )  6n/n  , 

D  o  o  Zn  2  o  on  op  o 

p  x„„  L 


L  — iSi 


»  • 


% 
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wit  h 


n  = 
s 


v  2  2 

E  Y 
o  *on 


2  2‘ 
16ira  mo 


(42) 


The  term  (1  -  iWn^)  6n/nQ  is  placed  on  the  right  side  of  the  equation  so  that 


an  iterative  solution  for  6n/nQ  may  be  obtained  by  using  the  Green's  function 


for  the  left  side.  Two  finite  temperature  solutions  for  6n  obtained  in  this  manner 


using  the  correct  finite  temperature  zero  order  profile,  nQ(r),  are  shown  in  Fig.  3. 


The  zero  temperature  solution  of  Eq.  (32)  is 


<5n 


"s[j5 


<xon  r/a)  "  (1  +  2kZ‘,2/*on)  Jo  (*on  r/a)  ]  cos2kz. 


(43) 


and  is  plotted  as  the  upper  curve  in  Fig.  3. 

On  axis,  the  ponderomotive  force  expels  particles  from  the  regions  where  |e 


is  large  and  causes  a  density  increase  where  |E  |  is  small  as  illustrated  in  Fig.  4. 

z 


The  electrostatic  field  arising  from  this  density  rearrangement  balances  the  pondero- 
motive  force.  Since  this  effect  is  two  dimensional,  the  charge  bunching  along  the 
axis  produces  a  fringing  electric  field  which  influences  the  force  balance  at  outer 


radii.  Since  6 | E^ |  decreases  monotonically  with  r,  a  radius  is  reached  at  which  5^ 


the  fringing  electric  field  due  to  charges  at  inner  radii  is  greater  than  the  field 
required  to  balance  the  ponderomotive  force  at  that  radial  position.  From  this  radius 
outward,  the  density  rearrangement  changes  sign,  i.e.,  density  increases  axially 


where  | E^ |  is  large  and  a  density  depletion  occurs  where  |Ez|  is  small.  This  effect 


is  seen  in  Fig.  3,  and  is  in  direct  contrast  to  the  situation  in  a  neutral  plasma  in 
which  this  fringing  effect  is  not  present;  in  neutral  plasmas  the  ponderomotive  force 
is  balanced  against  plasma  pressure  and  Poisson's  equation  is  never  utilized  because 
quasi-neutrality  is  satisfied. 

The  radial  dependence  of  the  nonlinear  density  modification  6n(z  =  0)  is 
shown  in  Fig.  3  for  zero  temperature,  K^a  =  25,  and  K^a  =10;  6n  has  the  same 


j 

•4 
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axial  dependence  as  6 | E  |  of  cos2kz.  At  z  =  0,  6 [ E  |  has  a  maximum  and  thus 

z  z 

fin  is  negative  on  axis  and  positive  at  outer  radii  as  shown  in  Fig.  3.  At  an 

axial  position  where  6 ] j  is  negative,  fin  becomes  positive  on  axis  and  negative 

at  outer  radii.  The  primary  effect  of  finite  plasma  temperature  occurs  in  the 

sheath  within  several  Debye  lengths  of  r  =  a.  For  zero  temperature  the  density 

drops  discontinuously  to  zero  at  r  =  a,  but  for  finite  temperature  the  density 

drops  smoothly  to  zero.  Finite  temperature  also  gives  rise  to  a  slight  reduction 

in  the  magnitude  of  fin  near  the  column  axis. 

It  can  be  concluded  that,  for  typical  nonneutral  plasma  parameters,  the  pon- 

2 

*  3  a  ~  2 

deromotive  force,  -  z  - - 6 1  E  |  ,  is  balanced  against  an  axial  electrostatic 

O  Z  /  Z  Z 

„  4rau> 

force,  -  z  —  qfi<j> ,  rather  than  the  much  smaller  pressure  force.  Pressure  forces 

only  play  a  significant  role  in  the  sheath  region  (the  outer  several  Debye  radii) . 

The  associated  density  modification  fin  is  obtained  from  5<J>  through  Poisson's 

equation;  the  force  balance  equation  is  of  little  use  because  of  the  near  exact 

cancellation  between  the  ponderomotive  force  and  the  electrostatic  force.  The 

modified  electrostatic  potential  6cf>  has  a  radial  electric  field  -  rr—  6<)>  asso- 

o  r 

dated  with  it  which  is  not  balanced  against  the  ponderomotive  force  since  the 
radial  component  of  the  ponderomotive  force  is  negligibly  small  in  the  strongly 
magnetized  limit.  This  radial  electric  field  changes  the  rotational  velocity 
of  the  plasma.  Substituting  Eq.  (40)  into  Eq.  (34)  gives  this  nonlinear  velocity 
modification  SV^ , 


6V0  =  Vs  VXon  r/a)  Jl(xon  r/a)  C°s2kZ’ 


„2_  2 

q  Eo 

4m2co2n* 


(45)  » 
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The  radial  dependence  of  SV  is  shown  in  Fig.  5.  A  shear  is  produced  in  the 
velocity  profile,  but  we  do  not  examine  any  possible  instabilities  resulting 


from  this  shear.  Since  both  6n  and  6V  have  an  axial  dependence  of  cos2kz, 
the  canonical  angular  momentum  is  on  the  average  unchanged,  as  expected  for  an 
axisymmetric  wave. 


Nonlinear  Frequency  Shift 


The  density  rearrangement  Sn  produced  by  the  ponderomotive  force  of  a  large 
amplitude  standing  wave  can  cause  a  nonlinear  modification  in  the  dispersion  rela¬ 
tion  of  the  modes  supported  by  the  plasma.  In  addition  to  the  self-modification 
of  the  large  amplitude  mode,  all  modes,  e.g.,  plasma  noise,  satisfying  proper  wave 
number  matching  conditions  are  nonlinearly  modified.  In  this  section  the  nonlinear 
resonant  frequency  shift  of  Gould-Trivelpiece  modes  resulting  from  externally  driving 
the  plasma  at  frequency  to  is  found.  For  simplicity,  all  calculations  in  this  sec¬ 
tion  are  performed  for  zero  temperature  and  an  initially  flat  density  profile,  since 
it  was  demonstrated  in  Sec.  Ill  that  thermal  effects  are  small  for  typical  nonneutral 
plasma  parameters. 

An  idealized  antenna  that  excites  waves  with  a  given  axial  wavelength  can  be 
obtained  by  cutting  the  conducting  cylinder  into  sections  of  length  X/2  and  driving 
adjacent  sections  with  alternating  polarities,  as  illustrated  in  Fig.  6.  Waves  of 
axial  wavenumber  k  =  2tt/X  and  odd  harmonics  thereof  are  excited  in  the  nonneutral 
plasma  column.  The  linear  plasma  response  is  given  by  the  superposition  of  modes 
which  satisfy  the  proper  boundary  condition  at  r  =  a. 


4V 


rf 


.L  TT  (2j+l) 
J=o 


j  [(2j+l)(u)  2/to2-l)1/2kr] 

- ^ -  sin[  (2j+l)kz ] , 

Jot(2j+l)(mpZ /(/ -l)1/Zka] 


(46) 


in  which  V  -  is  the  amplitude  of  the  applied  voltage.  In  order  to  make  the 


plasma  response  finite  when  the  excitation 

normal  mode  frequency,  collisional  damping 

2  2  2 
of  it)  /aiOo+iv)  -  1  for  a)  /w  -  1,  where  v 
P  P 

assumed  small  compared  to  u.  If  we  assume 


frequency  is  in  the  vicinity  of  a 
is  included  through  the  substitution 
is  a  collision  frequency  that  is 
that  <j o  is  near  a  resonance  of  the 


fundamental  k  mode,  i.e.. 


then  we  mav  separate  the  peak  amplitude  $  of  this  mode  from  Eq.  (46)  as 


*  X  J,  (\)  [(to  -  to  )  +  v  /4] 

on  l  on  o 

in  which  is  defined  through 

4>  =  4>  J  (x„„  r/a)  sinkz. 
p  o  on 

This  large  amplitude  wave  produces  a  nonlinear  density  modification  given  by 


Eq.  (43): 


6n/nQ  =  e(r)  cos2kz, 


e(r)  = 


,2  2  2 
k  <p  x 
p  on 


•  2  2 

16-rrn  min  a 
o  o 


i_  [j2 
J1 


(XQn  r/a)  -  (1  +  2k2a2)j2(xon  r/a) 


Next  we  calculate  the  nonlinear  modification  to  the  Gould-Trivelpiece  dispersion 
relation  for  small  amplitude  test  waves  supported  by  this  modified  plasma  equi¬ 
librium.  The  partial  differential  equation  for  $ ,  Eq.  (15),  when  expressed  in 
the  strongly  magnetized  limit  and  with  the  time  derivatives  retained  is 


„2  ,  ,  3  23.  _ 

1  V  *  +  dz  “p  Tz  *  =0, 

-  V.  1  ✓ 


The  fractional  density  change  6n/n  is  included  in  Eq.  (52)  by  substituting 

2  2 
to  (1  +  Sn/n  )  for  to  ,  giving 
p  o  P 


32  2  A  ^  2  J‘ 

- ^  V  <J>  +  to  — 


=  -wd2  e(r)  h  fcos  2kz  b 


2 


(53) 


Equation  (53)  is  solved  for  the  eigenmodes  of  the  modified  system  when 
6n/n  «  1  by  a  perturbation  technique  in  which  the  solutions  are  expanded  in 
terms  of  the  solutions  of  the  equation  with  e  =  0.  Assuming  that  the  solution 
consists  of  a  perturbation  to  a  small  amplitude  test  wave  of  the  form 


=  f(z,t)  J£(x£m  r/a)  exp (i£0 ) , 


in  which 


f(z,t)  =  <|>  exp{i(k.z  -  o).t  -  6m.  t)  ]  +  exp[i(-k.z  -  to  t  -  6co.t)],  (5! 

Ill  111 

and  kj  and  to^  are  related  through  the  dispersion  relation  for  the  unperturbed 

system,  2 

,22  X£m 

y  -  -  - — r—2  • 

1  -  (0  /<0  . 
p  1 

The  e(r)  term  causes  the  frequency  shift  6ok  appearing  in  Eq.  (55)  and  introduces 
additional  radial  and  axial  modes  with  frequencies  other  than  w.  In  this  study  we 
do  not  consider  these  additional  modes;  instead  we  concentrate  on  the  frequency 


shift  6m j . 

Inserting  Eqs.  (54)  and  (55)  into  Eq.  (53),  multiplying  by  rJa(XJm  r/a)exp(-i£0) , 
integrating  in  r  from  0  to  a,  and  projecting  out  the  z  independent  component  gives 


(6uj/«j)  a2  J£+1(xim)[l<f>+|2  +  k_|2]  =  -($!!  <fr+  +  ^_^*)jrdr  e(r>  j2(Xj 


r/a).  (57) 


The  e(r)  term  introduces  a  coupling  between  the  right  going  and  left  going 
waves  when  k^  =  k.  The  physical  reason  is  as  follows:  The  density  ripple  has 
an  axial  wavenumber  of  2k  while  the  difference  in  wavenumber  between  left  traveling 


and  right  traveling  waves  is  2k^ .  Thus  when  k  =  k^  the  density  ripple  couples 
the  left  traveling  and  right  traveling  waves. 
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Specializing  to  an  axially  standing  wave  with  =  $./2i,  $ 


-<j> .  /2i ,  and  — 

J  » 


k  =  k  ,  Eq.  (57)  becomes 


<5w  .  /w  .  = 
J  J 


r  dr  £(r)  J£(x£m  r/a>)(a2  J£+l  (x£m})  * 

o  '  ' 


(58) 


and  <j> .  is  a  constant  amplitude  of  the  small  amplitude  test  wave, 
into  Eq.  (58)  gives 

-  k2|U  2  2  2  2  2 

=  - w - -  Ux„n  +  2k  V)  C  -  X„n  D], 

J  J  16ira2n  mmV  (x#  ) 
o  o  £+1  A£m 


Substituting  Eq.  (51) 


where 


I  *i 


udu  J2(Xn  u)  J2 (x  u), 
0  £  A£m  o  Aon  ’ 


D  s  judu  J*(xtm  u)  j2(xon  u) 


Since  Eq.  (59)  describes  the  nonlinear  frequency  shift  of  axially  standing 
small  amplitude  test  waves  (of  linear  frequency  w  )  due  to  the  presence  of  the  large  1 

yy. '  \  -  '■] 

amplitude  axisymmetric  standing  wave  at  frequency  we  require  that  either  m  ^  n  ^  ^ 

or  £  ^  0.  This  frequency  shift  is  nonzero  only  when  the  test  wave  has  the  same  axial  1 


dependence  of  sinkz  as  the  large  amplitude  wave. 

The  previous  analysis  also  describes  the  self-modification  of  the  large  amplitude 


mode.  If  m  =  n  and  £  =  0  then  w .  becomes  u>  ,  6u> .  becomes  6w  ,  d> .  becomes  <$>  ,  and 

J  o’  j  o’  Yj  p’ 

Eq.  (59)  describes  the  nonlinear  frequency  shift  of  the  large  amplitude  mode: 


(Soj  /u>  - - 

°  °  18ira2 


-*2IU2 


■— _ rx2  e  + 

n  mu)2 J2  (X  )  L 
o  o  1  on 


2  2 

k  a  F 


where 
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E  = 


F  3 


udu  ^(XGn  u)  3-Xon  u), 


”du  Jo(xon  uK 


(62) 


Replacing  u)o  in  the  resonant  denominator  of  Eq.  (48)  with  +  ^WG»  usin8  Eq.  (61) 
for  5ojo,  and  defining  the  scaling  parameters. 


8oi  V  c 
o  rf 


™XonJl(xon} 
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8k  a)  V  c 

_ o  rf _ 

3  2  3  2  2 

*  3  nomV  XoJl(Xon) 


2  2  2 
'XonE  +  ka 


(63) 


(64) 


Aw  =  a)  -  a)  , 
o 


(65) 


gives  a  nonlinear  equation  for  the  wave  amplitude, 


iy*.i 


f2Au)/v  +  2p|(J>  / d>  |^]^  +  1 
P  s 


(66) 


The  resulting  scaled  wave  amplitude  c(> ^ /<J> ^  obtained  by  numerical  solution  of 
Eq.  (66)  is  plotted  in  Fig.  7  for  several  values  of  the  nonlinearity  parameter  p. 
For  values  of  p  >  1  the  plots  of  wave  amplitude  become  multiple-valued.  Since  only 
one  value  is  physically  possible,  this  causes  the  system  to  follow  different  ampli¬ 
tude  versus  frequency  curves  depending  upon  the  direction  the  excitation  frequency 
is  swept  through  the  resonance  frequency.  Figure  8  shows  this  hysteresis  effect 
for  p  =  2.  If  the  frequency  w  is  swept  slowly  up  (i.e.,  increased)  through  reso¬ 
nance  compared  to  the  system  relaxation  time  1/v,  the  system  response  follows  the 
curve  from  A  to  C,  jumps  up  to  point  E,  and  continues  along  the  curve  to  F.  When 
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the  frequency  is  swept  down  slowly,  the  system  response  follows  the  path  FDBA. 
An  analogous  hysteresis  effect  has  recently  been  experimentally  observed7  in  a 
two-dimensional  nonneutral  plasma  consisting  of  charges  suspended  below  the 
surface  of  liquid  helium. 
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V.  Conclusions 

The  present  analytical  study  illustrates  the  principal  features  of  the 
ponder omotive  effect  in  nonneutral  plasmas.  Whereas  in  a  neutral  plasma  the 
ponderomotive  force  is  balanced  against  plasma  pressure,  in  a  nonneutral  plasma 
it  is  balanced  against  self-consistent  space  charge  fields.  The  consequence  of 
self-consistency  in  a  finite  geometry  is  that  the  local  density  change  produced 
by  the  ponderomotive  force  is  not  simply  related  to  the  local  strength  of  the 
electric  field  amplitude,  as  is  the  case  for  a  neutral  plasma.  Consequently, 

2 

the  usual  formula  6n/n  =  —  |  E I  /16Trn  T,  does  not  hold  in  the  nonneutral  system. 

o  _  o 

Instead,  the  global  properties  of  the  system  are  important.  This  leads  to  a 

situation  in  which  near  the  axis  of  the  nonneutral  column  a  density  change 

given  approximately  by  6n/n  =  -(y  /a  )  6 1 E  I  /8irn  mw  ,  in  which  -  x  / a  is 
J  J  o  on  1  z1  o  on 

-  2  ( 

the  radial  curvature  of  I E  I  on  the  axis.  At  outer  radii  the  radial  curvature 

1  z 1 

-  2 

of  I E  |  becomes  positive,  and  a  reversal  in  sign  of  the  density  change  results. 
This  effect  is  associated  with  the  fringing  of  the  space  charge  field,  as  dis¬ 
cussed  in  Sec.  III. 

In  addition  to  the  density  changes,  the  ponderomotive  force  is  found  to 
alter  the  zero  order  rotation  of  the  nonneutral  plasma  column.  In  particular 
it  introduces  a  velocity  shear  whose  consequences  have  not  been  explored  in 
this  study. 

To  obtain  concrete  results  which  are  manageable  analytically,  the  present 

f 

study  has  examined  in  detail  the  ponderomotive  changes  produced  by  an  ideal 
large  amplitude  standing  Gould-Tr ivelpiece  mode.  For  simplicity,  the  plasma 
is  taken  to  be  infinite  in  axial  extent.  While  this  model  is  necessarily  re¬ 
stricted,  the  results  obtained  are  expected  to  be  characteristic  of  these  systems. 

/ 

However,  the  usage  of  more  realistic  waveforms  and  consideration  of  finite  axial 
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extent  should  yield  more  accurate  results  that  may  fit  closer  the  actual  experi¬ 
mental  situation.  The  additional  anticipated  effect  not  included  in  the  present 
model  is  an  axial  elongation  of  a  finite  length  column. 

This  study  has  also  evaluated  the  lowest  order  consequences  of  the  nonlinear 
changes  in  density  produced  by  the  ponderomotive  force.  It  is  found  that  the 
linear  dispersion  relation  is  modified,  and  results  in  a  frequency  shift  for  some 
of  the  standing  waves  in  the  system  which  satisfy  the  appropriate  matching  condi¬ 
tions  explained  in  Sec.  TV.  The  result  is  evaluated  for  additional  test  waves  as 
well  as  for  the  large  amplitude  wave  responsible  for  the  density  rearrangement. 

In  this  latter  case,  an  interesting  hysteris  effect  is  obtained.  It  arises  because 
the  resonant  response  function  (i.e. ,  the  usual  Lorentzian  line)  becomes  amplitude 
dependent.  Consequently,  when  the  nonlinear  frequency  shift  becomes  comparable 
to  the  collisional  broadening  of  the  resonance,  an  asymmetry  in  the  response  re¬ 
sults.  As  the  external  driving  frequency  is  swept  slowly,  different  amplitudes 
can  be  generated  which  correspond  to  the  various  solutions  of  a  cubic  equation. 

An  analogous  effect  has  been  recently  observed  in  an  experiment ^  involving  a  two 
dimensional  nonneutral  plasma  generated  by  attaching  charges  below  the  surface 
of  liquid  Helium.  The  present  results  indicate  that  a  similar  effect  should  also 
occur  in  a  3-dimensional  nonneutral  plasma. 

Finally,  the  feasibility  of  experimental  observation  of  the  effects  predicted 

is  assessed.  Considering  a  presently  attainable  peak  electron  density"* 

7—3  8 

n  =  1.3  x  10  cm  (w  =  2  x  10  rad/sec)  and  column  radius  a  =  2  cm,  one  finds 
P  P 

a  relation  between  the  peak  potential  amplitude  $  for  the  standing  wave,  taken 

as  the  lowest  radial  Gould-Trivelpiece  mode,  and  the  fractional  density  dn/n^, 

1/2 

<t>  -  45(<5n/n  )  volts.  The  self-consistent  A  can  be  related  to  the  external 

P  o'  p 


For  a  Q  -  100,  an 


potential  V  through  the  Q  of  the  mode,  i.e.,  <J>  =  2QV  . 

Kr  p  Kr 

external  excitation  of  the  order  of  22  mV  is  predicted  to  result  in  a  1%  density 

change,  which  should  be  measurable.  Concerning  the  distortion  of  the  resonant 

curve  due  the  nonlinear  frequency  shift,  it  should  be  noted  that  interesting 

effects  begin  to  appear  when  the  nonlinear  parameter  p,  defined  in  Eq.  (65),  is 

of  order  unity.  The  condition  p  =  1  requires  an  external  excitation 
-3/2 

V  =  70  Q  volts  which  again,  for  a  Q  =  100,  predicts  Vnt,  -  70  mV,  thus 
Kr  Kr 

suggesting  that  the  effect  should  be  amenable  to  experimental  investigation. 
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FIGURE  CAPTIONS 


a)  Schematic  of  geometry  and  large  amplitude  £  =  0  axially  standing 
wave,  b)  Typical  zero  order  density  profile  for  a  =  10  and 
radial  dependence  of  the  lowest  order  Gould-Trivelpiece  mode. 

Radial  dependence  of  6 ] |  for  the  lowest  order  radial  mode,  evaluated 
at  the  axial  position  z  =  0.  The  axial  dependence  is  cos2kz. 

Radial  dependence  of  the  nonlinear  density  rearrangement  Sn/n^, 
due  to  the  mode  in  Fig.  2,  evaluated  at  the  axial  position  z  =  0. 

The  axial  dependence  is  cos2kz.  6n  is  shown  for  three  different 
plasma  temperatures  in  the  limit  ka  <<  1. 

Schematic  of  the  plasma  response  to  the  ponderomotive  force  of  a 
standing  wave  showing  the  fringing  electrostatic  field  lines  for  a 
positive  ion  plasma.  The  physical  interpretation  follows  Eq.  (43). 
Radial  dependence  of  6V  evaluated  at  the  axial  position  z  =  0. 

The  axial  dependence  is  cos2kz. 

An  idealized  antenna  for  exciting  standing  waves. 


Scaled  wave  amplitude  (ji  /<j>  versus  excitation  frequency  w  for  three 

P  s 

values  of  the  nonlinear  parameter  p. 

Scaled  wave  amplitude  $  /(j>  versus  excitation  frequency  io  for  p  =  2. 

P  s 

When  frequency  is  swept  up  (i.e.,  increased)  the  wave  amplitude  follows 


ACEF.  When  frequency  is  swept  down  (i.e.,  decreased)  the  wave  amplitude 


follows  FDBA. 
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